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Abstract 

The CDW Hartree-Fock state at half filling and half electron per unit cell is 
examined. Firstly, an exact solution in terms of Bloch-like states is presented. 
Using this solution we discuss the dynamics near half filling and show the 
mass to diverge logarithmically as this filling is approached. We also show 
how a uniform density state may be constructed from a linear combination 
of two degenerate solutions. Finally we show the second order correction to 
the energy to be an order of magnitude larger than that for competing CDW 
solutions with one electron per unit cell. 
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I. INTRODUCTION 

For almost two decades the fractional Quantum Hall effect has been one of the main focus 
of interest in condensed matter physicsua. Much progress has been made towards its un- 
derstanding, and even general theories exist today aspiring to be fully coherent descriptions 
of the underlying physics!. These theories do not rest on an actual solution of the basic 
quantum mechanical equations of motion, however, but are rather cast from ansatz wave 
functions exhibiting a large overlap with accurate numerical solutions of such equations for 
a few particles. A more basic theory, although highly desirable, is very difficult to attain 
because electron-electron interactions and correlations are at the core of the effect. 

A first step in a perturbative approach was developed in the early years by Yoshioka and 
Lee, who constructed a mean field Hartree-Fock theory for the spin polarized case (HFTl)i. 
It received little attention, however, because it failed to provide the empirical selection 
rule that distinguishes even from odd denominator filling fractions, which characterizes the 
effect in an essential way. It was later shown by one of us, that if the mean field solution 
is constructed in a slightly different way, such distinction arises naturally, giving the same 
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gap structure at the various filling fractions as experiment, and producing the proper step 
like Hall conductivity dependence with magnetic field (HFT2)i'i. It predicts a gap at every 
odd denominator fraction, and a metallic state at all even denominators, within a band 
spectrum whose fine structure at or near the Fermi energy scales with the denominator 
in the fraction. Energies were too high, however, even compared to the state proposed 
by Yoshioka and Lee. One potential problem of Hartree-Fock spin polarized states in the 
lowest Landau level is that they necessarily have non uniform electron density at fractional 
filling!, and, if space fluctuations are severe, possibly pin the many-body electronic state to 
the underlying impurities, an effect for which there appears to be no experimental evidence 
except perhaps at fillings below 1/7. However, it should be also stressed that a crystalline 
ordered state ("Hall Crystal") can be fully compatible with a quantized Hall conductivity 
thanks to its magnetic field dependent crystal parameters^. 

Simple charge density wave (CDW) mean field states are found assuming they form a 
periodic lattice of rectangular or triangular geometry. The size of the unit cell is an additional 
degree of freedom characterized by the quantity 7, the number of electrons per unit cell. The 
main difference between the mean field theories described above is that HFT1 assumes one 
electron per unit cell, that is, 7 = 1, while HFT2 assumes 7 = 1/2 or some other fraction of 
denominator 2 k , with k an integer. Keeping the geometry fixed and changing 7 yields local 
minima at 7 = 1/2, 1. Although neither show energies near the one corresponding to the true 
ground state, one may ask which one is a better perturbative precursor of the ground state. 
HFT1 gives a lower energy, but HFT2 provides the basic selection rule distinguishing even 
and odd denominators as required by experiment, and has less pronounced charge density 
fluctuations. 

In this work we present an analytic solution for the v = 1/2, 7 = 1/2 state. In particular, 
we address the question of sensitivity to perturbation theory in both theories. Within HFT1, 

2 

Yoshioka and Lee obtained a correction to second order of about 0.002 in units of — , where 
r a is the magnetic length@. Here we show the same correction to be an order of magnitude 
larger for HFT2. This result confirms the suggestion made in a former work by one of 
us about that this latter state is more sensitive to correlations than HFTli. Our solution 
is constructed in terms of Bloch-like running waves which solve the Hartree Fock problem 
exactly^ and form a complete orthonormal set, save for a single point in the Brillouin zone0. 
This remarkable property may be unique to filling 1/2 and 7 = 1/2, since then exactly one 
flux quanta traverses the unit cell of the underlying CDW . Having analytic solutions allows 
for a study of the dynamics at the Fermi surface. We find that the cyclotron effective mass 
diverges logarithmically as half filling is approached, in agreement with RPA estimates^ 
and with some experimental dataEl. 

In Section 2 the single particle solutions and their self-energies are given. Close expres- 
sions in terms of products of elliptic theta functions are presented, and the order parameter 
is explicitly given. We also show that degenerate Hartree-Fock solutions may be superposed 
to form a state of uniform charge density. In Sec. 3 a semiclassical analysis of the response 
to a small additional field is given, and the resulting cyclotron mass discussed. Section 4 
is devoted to the evaluation of the second order correction to the energy. A short review 
and discussion of our results is given in the Summary. In the appendices we derive some 
symmetry properties of our solutions and show that filling the rotated central square of the 
Brillouin zone indeed yields the lowest self-consistent energy solution. 



2 



II. ANALYTIC SINGLE PARTICLE SOLUTIONS AND ORDER PARAMETERS 



We consider N electrons constrained to move in a plane of area S under a perpendicular 
magnetic field B. As shown by Wannier0, to study this system one can construct an orthog- 
onal set of single particle Bloch-like states from the zero angular momentum eigenfunction 
in the lowest Landau level 
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Eigenstates are labeled by a wave vector k = p/h where p is the quasi- momentum, and the 
sum runs over all integers £%, £2 defining a planar square lattice L, with £ = a(i\, £2, 0), a = 
2Ttr 2 and r Q = y^^.The magnetic translation operator acting on any function / intro- 
duces a phase: 

T, /(x)= exp(— A(/).x) /(x - £), (3) 

where the vector potential is assumed in the axial gauge A(x) = ^(— X2, xi, 0) and the 
electron charge e is taken with its negative sign. 

We want to study the case in which one flux quanta traverses each lattice cell. Then, 
there is one state per plaquette in the lowest Landau level and at filling one half the charge 
in each cell is just half the electron charge, and thus, 7 = 1/2. The wave-function (2) and the 
normalizing factor can then be expressed in terms of elliptic theta functions thanks to 
their simple properties under special shifts of the complex arguments in their quasi-periods. 
They are given, respectively, by 
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where r = | . The components of the wave vectors appearing in these expressions are 
related through k* =k+(nxx-i x)/2rg, where n = (0, 0, 1) is the unit vector normal to 
the plane. The orthonormal set is well defined except at the single point k = (ir/a, ir/a, 0), 
where the norm vanishes^!. In what follows we shall ignore this singular point. 
As may be easily verified, the functions y?k satisfy the eigenvalue equation 

T t p k (x) = X k {£) <^ k (x), (4) 
\ k (£) = (-l) £l£a exp(-z k.£). 

As was shown earlier they are exact solutions of the Hartree-Fock (HF) single particle 
Schrodinger equation associated to an arbitrary Slater determinant formed by selecting an 
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also arbitrary group as filled stateJS. This occurs because the HF single particle Hamil- 
tonian commutes with all translations leaving L invariant!. Since the functions (2) are 
common eigenfunctions of the commuting magnetic translations leaving invariant the lattice 
L and the set of eigenvalues (4) uniquely determines them, the HF potential associated 
with the Slater determinant can not change those eigenvalues. Therefore tp^ should be an 
eigenfunction. 

The explicit expression for the self-energy of the state cp^ has the form,! 

<W - ^I>(Q) A-(Q,e XP (-^) (H£W-¥> - k 

(5) 

where Q = 27r(n 1 , n 2 , 0)/a is the set of reciprocal lattice vectors and I is the modified 
Bessel function. The momentum dependence of the self-energy is fully contained in the 
factor Ak(Q) = 27rr 2 Ak(x*)/S', where x* = r 2 Q n x Q. The order parameter is in turn given 
byi, 

27rr 2 

A(Q) = E A k (*) (6) 

D kGF 

where the sum is over all filled states in the Brillouin zone B, the set of which we call F C B . 
From general symmetry properties it is shown in Appendix B that an energy minimum is 
obtained among all possible Slater determinants by filling states inside the square bounded 
by the constant energy lines ±k x ± k y = n/a, which we take as the Fermi surface^!. The 
number of states in this square is just half the total, as required, since we are studying the 
half filling case. Note that the energy function is continuous across these lines so there is no 
gap in the single particle spectrum at the Fermi energy. Turning the sum into an integral 
one then gets, 



A(Q) = (-lf^a 2 / y^expH k.x*) (7) 




_ _-_ y sin[(wi - 7i 2 )f] sin[(wi + n 2 )f] \ 
7r 2 e^o n\ — n 2 , + e J 

Notice that these quantities vanish if n\,n2 have the same parity, save at the origin. Using 
the above expressions the total energy per electron may be easily computed, to obtain 
ei = —0.394—, a value somewhat above the result for 7 = 1, ei = — 0.443— ? . 

The charge density may also be obtained from the above expressions since their Fourier 

r 2 Q 2 

coefficients and (7) are related through A(Q) = 27rr 2 e _2 ^p(Q), yielding a space fluctuation 
of a mere 17% about the average no = N/S, while the 7 = 1 state electron density varies 
between 0.1n and 1.85n , representing a fluctuation of close to 90% about the average.0 

Finally in this section we would like to underline the curious point that using our Hartree- 
Fock solutions it is possible to construct a state of uniform charge density. To see this we 
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consider the function \l/ = (Q HF + Q HF ) / y/2, where as before <& HF is the Slater determinant 
formed with all states in F, while Q HF is the Slater determinant of all states in B — F, 
the complement of F in the Brillouin zone B. Both regions are separated by the Fermi 
surface, the square ±k x a ± k y a = tt. Because all single-particle states are orthogonal to 
each other the crossed term in the charge density vanishes and one has p(x) = (p HF (x) + 
p HF (x))/2 = J2k&B |^k(x)| 2 /2, half the density of a completely filled Landau level, which 
equals (47rr2) _1 and is uniform throughout the sample. 

We now show that the above linear combination has the minimal energy. To see this 
we first note that § HF is obtained simply by shifting all momenta in § HF by the vector 
5 = (1, — l)7r/a, since in the extended zone scheme this vector maps all states in F onto 
B — F . However, according to Eq. A. 2 such shift may be done with the aid of a single 
magnetic translation Tr,, with R = r 2 n x 5 =(1, l)a/2. Thus, both Slater determinants 
are related by a single translation in space. But, all the N states forming each of these 
Slater determinants are orthogonal with all the ones defining the other, and also the energy 
operator is a linear combination of products of merely four creation or annihilation operators. 
It therefore follows that the mean value of the energy operator in the state \I/ is just half 
what one obtains adding the mean values of & HF and & HF , and thus the mean energy of 
the superposition state coincides with the HF energy. 



III. DYNAMICS AND MASS 

It has been suggested that the effective mass of the charged carriers may diverge at the 
Fermi surface, a property that remains controversial. Taking advantage of our analytic 
results we examine these questions. The semiclassical equations of motion may be solved for 
particles at the Fermi surface moving in the external field corresponding to half filling, Bi, 
with the result that if motion starts at k = ^(1,1) then, after a time t, the particle quasi- 
momentum has reached the point k x = - arctan[exp(— t/T)], ak y = 7r — ak x , with a time 
constant of order r = 2h/e*, where e* = — . It thus takes an infinite time to reach a corner 
of the square Fermi surface, and the trajectory in real space is a straight line covering just 
a fraction of the unit cell. One can thus claim that the Fermi particles behave as if there 
was no external field at all, in agreement with previous workS 

In order to examine the effective mass that this solution provides we consider the pe- 
riod of a cyclotron orbit when the filling is slightly above or below 1/2. For definiteness 
we assume the later. For simplicity me keep the lowest order Fourier components in the 
dispersion relation only, and assume that the slight change in filling fraction does not 
perturb significantly the particle density. The self energy of state k is then of the form 
e = e Q — ei(cos(a/c a; )+cos(a/c 3/ )), with e Q = — \J\ e and e x = 0.087 e*. We study the dynamics 
governed by the equations 

.dk e de „ 

n — = — — x B. 

dt he ok 

For an orbit over the Fermi surface, no longer square in shape, this equation becomes, for 
the component k y , 

dk v ae\ 
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from whose solution and the constant energy condition the itinerary of k x may be extracted. 
Here 5 2 = Tic/eB = 2-nr 2 , and 77 = (e — e )/e\ measures the departure from the half filling 
square Fermi surface. Integrating this equation one obtains for the time it takes for a particle 
to go around the energy contour once, 



TTCi 77 

where K(u) denotes the complete elliptic function of the first kind. A cyclotron effective 
mass may be obtained from this expression through the usual definition m* = eTB /(2irc). 
One obtains the result m* = 8h 2 n K ( n -^) / (ne^) . Noting that K( n -^) ~ — 0.57/7 l°g V as 
77 — > we get the final result, 

16.7 h 2 n 

m = log (ry), 

e* 

which diverges logarithmically as half filling is approached. 



IV. ENERGY PER PARTICLE IN SECOND ORDER 

We now turn our attention to the second order correction to the energy per particle 
for the v = |, 7 = ~, state. A similar evaluation for 7 = 1 was done by Yoshioka and 
Lee, obtaining; a correction of the order of 0.5% of the total result.! . However, as it was 
underlined inB, the increased degree of overlap of the single electron states associated to the 
7 = | wave-functions could change the situation drastically. The study of this question is 
the main objective of the present paper. 

For evaluating the second order correction the following formula will be employedB 

r(HF)\ I rh A 1 I (U _ ui HF )\ I *HF\ 



E m = I (H - Hf F >) I ^) E(HF) _ E S ^ I (H - Hf F) ) I S> HF ) (8) 

I (§HF \H\&i) | 2 



E 



E( HF ) - Ei 



where Q HF , E^ HF ^ are the Slater determinant and total Hartree-Fock energy, respectively, 
and H is the projection of the exact Hamiltonian onto the first Landau level. The many 
particle excited states $j are Slater determinants constructed with the Hartree-Fock basis 
states {<£>k} , k G B. It follows that (Q HF | $j) = 0, a property that allowed to write 
the last equality in (8). In the second quantized representation the Hamiltonian H will 
have non-vanishing matrix elements linking the HF state and excited states of the form 
I $j) = a v a v ' a^, \ & HF ), where a|" creates an electron of wavevector £, etc. The index i is 
a shorthand notation for pairs of filled (77, 7/ G F) and empty (£, £' G B — F) electron states. 
The total energies of the excited states are given by E{ = E^ HF ^ + e(£) + e(£') — £(77) — e(rj'). 
Then, the second order correction can be rewritten in the form 

F (2) = V V \^ HF \ Ha V a rf44\® HF ) I" (a) 
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where the total projected Hamiltonian has the form 

H= — j f rfxrfx' #*(x)**(x') | 1 —— *(x')*(x) (10) 

2 J J 



e 2 

tEE M(a, a' | (3, (3') a+ a+, , 

Z a,a> /3,/3> 



with the matrix element of the coulomb interaction given by 

M(a,a' | (3,(3') = J J rfxrfx' ^ ]_ ^ ^(x^x). (11) 

By using the anti-commutation relations [a q , aj] = <5 q , q ', formula (11) can be expressed in 
terms of our basis as, 

F(2) = v v |e 2 //^x^; Y (x,xQ^^(x,xO | 2 

^)+ C (^)- C (o- C (eo ' 1 j 

where the two particle states <3> are defined as 

(pJx)(p ,(x') — ld ,(x)o9 r) (x') 
$W(x,x') = ™ ^ ; ^ ^ ; . (13) 

This quantity can be evaluated by use of the series (2), to obtain 

M(a',a\ (3,(3') = N 4>0 5 { p\a + a', (3 + (3') I((3' - a, (3' - a') * (14) 
<(-r 2 n x (^-q))^(r 2 n x (/?' - a)) 

where the Coulomb interaction is contained in the function /, given by 



Q 



/?' - a + Q 



exp - a + Q) 2 ) exp (i r 2 n x (a' - (3') • Q) . 



(15) 



The inherent translation invariance of the problem is reflected in the delta function express- 
ing the conservation of the quasimomentum of the four particle states defining the matrix 
element, 

<$°fa W, £ + a = E $ (K) (v + + + Q). (16) 

Q 

Finally, by using the above expression for the matrix element, the second order correction 
to the energy per particle can be expressed as 



1 r r r dq dq' dp M(q, q' | p) 



e 



2 



N 2(2vr) 2 J J J (2tt) 2 (2vr) 2 (2tt) 2 e(q) + e(q') - e(p) - e(q + q'-p)er ' ( ' ' 1 
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where the momentum variables have been rescaled to be dimensionless through the changes 
of variables q = a 77, q' = a 7/, p = a £, and the usual equivalence in the high area limit 
SqS'(p) = S J j0^g{v) has been employed. The expression for the kernel in (13) becomes 

> w i i \ i w / ,„ , nx (q'-p). „ ,nx (q'-p), 
M(q,q'|p)=|/(q'-p,q-p)tf q ( ^ 



2tt y " q v 2tt 



r / h, / nx(p-q) nx(p-q) 2 1 

I(p-q,p-q)^q( = )#q'( o n / n xn / n xn n TTyT ( 18 ) 

2tt 2vr tf q (0)tf q /(0)tfp(0)tf q+q >_p(0) 



Relation (19) allows an estimate of e^ 2 ' by evaluating the three momentum integrals. The 
integration was performed partitioning the Brillouin zone into a lattice of (2n + l) 2 points, 
over which the integration variables take the values q(m 1; rri2) = 2^+1 ( m i ~ m 2 ; wii + ^2), 
— n < m 1 ,m2 < n, and similarly for q'(m 1 ,m2), p(mi,m 2 ). These partitions have the 
property that the points do not touch the boundary of F, and from which they are at least 
a distance equal to half the distance between the points of the partition. This property was 
implemented in order to avoid in a regular manner the singularity which appears when the 
momenta of the two states inside and outside F are all on the Fermi boundary. When at 
least one of the states is outside F the difference of energies in (53) is always non vanishing 

2 

due to the form of the self energy surface. The results in units of j^- for partitions involving 
up to 225 points are 
n e® (g) 

1 -0.0572742 

2 -0.02460706 

3 -0.02515106 

4 -0.02762403 

5 -0.02992157 

6 -0.03175841 

7 -0.0307658 

Comparison with the 7 = 1 value 0.002e 2 /£:r o shows an increase in the correction for 
7 = 1/2 of over an order of magnitude. Together with the significantly smaller space 
fluctuations of the charge density, we interpret this result as a larger sensitivity of the 
7 = 1/2 state to the introduction of correlations, leading to a faster lowering of the energy 
and possibly melting of the charge density wave. Extension of this work to other filling 
fractions will be reported elsewhere. 

V. SUMMARY 

An analytic solution of the Hartree-Fock problem at filling v = \ and half a particle per 
unit cell (7=5) has been discussed. The same state was formerly studied numerically.^ 
The solution is found to have a more uniform charge density in space than the 7 = 1 state 
of Yoshioka and Lee, and the correction to the energy is an order of magnitude larger than 
that obtained with one electron per cell. Besides yielding no gap as required for filling 1/2, 
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our results suggest that the Hartree-Fock state with half of an electron per unit cell is a 
better perturbative precursor to the true ground state of the system, and calls for a more 
detailed investigation of its properties for the filling considered, as well as other fractions. 

It should be stressed that after finishing this work Dr. N. Maeda have communicated us 
about two references:0'0 in which the same HF electronic state at v = | was considered by 
him and other collaborators as possibly related with the composite fermions at this filling 
fraction. In the present work, from a technical point of view, we have intended to present and 
analytical solution to the mean field problem as a new result. Alternatively, on the physical 
side, our main purpose was to decide about wether the enhanced electron overlapping at 
7 = | is able to also increase the correlation energy of the state at 7 = | over the one for 
7 = 1. 
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VII. APENDIX A 

The following properties of the basis functions under translations and reflections are 
needed in the body of the text. 

A. Translations 

Let us argue that the effect of an arbitrary translation on the basis functions is, modulo a 
phase factor, equivalent to a shift in the momentum label0. Operating over a basis function 
twice with the translation operator involving a lattice vector t and an arbitrary vector a one 
gets, 

T a T £ ^ p (x) = A p (£)T a ^ p (x) (A.l) 
= exp(-2— A(a)i) T £ T a y? p (x), 

where we have used Eq. (4) and the identity 

T a T e = exp(2— A(a)i)T,T a . 
he 

Again using (4), 

T e T a ^ p (x) = A p+2 « A(a) (£)T a y? p (x). 
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Then, taking into account that the set of eigenvalues defines uniquely the wavefunctions 
modulo a phase, it follows 

T a v? P (x) = F p (a) W 2 ^ A (a)(x). ( A - 2 ) 

That is, a magnetic translation is equivalent to a shift in the quasi-momentum. The factor 
F, being a phase, satisfies 

F p (a)F p *(a) = I. 

We now consider a special translation in the vector a = r 2 Q n x 5, with 5 = (—-, -)■ Using 
Eq. A. 2 one finds that it amounts to a shift in the momentum by the quantity 2^A(a) = 5. 
Then, the magnetic translation in a vector a transforms any momentum p G F into a 
corresponding momentum p 5 G B — F through the relation 

P 5 = p + 5 (A.3) 

thus transforming all of F in all B — F. 

B. Reflections 

Consider the axis D in the plane determined by the vector 

d=(l,l), (A.4) 
and define the operation R of a reflection on D within the plane by 

= (x 2 ,x 1 ) = R x, (A. 5) 

= R p, 

and its self- inverse property 

R = R- 1 . 

After acting with the operator R over any function / the following property follows: 

RT e /(x) = i?exp(-A(£).x) /(x - i) 

IP 

= exp(— A(£).Rk) f(Rx - I) 

IP 

= exp(— A(i^).x) /(i?(x - M)) 

he 



= T m Rf(x), 



from which one gets 



^p(x) = -L E ("I)' 1 ' 2 exp(^ i2p.€) T, 0(x) (A.6) 

meaning that the action of a space reflection on the basis function is equivalent to a reflection 
of its quasi-momentum. 



10 



VIII. APENDIX B 



That the lowest energy state for half filling is obtained by occupying all states within the 
square bounded by the lines ak x + ak y = n may be shown as follows. First, let us inspect the 
commutation properties of the reflection R under the line containing the vector d = (1, 1), 
as well as the particular translation T a which transforms the filled states F into the empty 
sates in B — F, with the Hartree-Fock single particle Hamiltonian. Both operations are 
defined in Appendix A. As the action of this hamiltonian is represented by a kernel with 
two arguments due to the projection on the first Landau level, it turns out to be useful to 
employ, when appropriate, the Dirac notation. Then the single particle HF equation can be 
expressed as 

H^ F) | y? P ) = e(p) | tp p ) (B.l) 
where the kernel associated to the Hartree-Fock Hamiltonian is 

2 

< F) (x,x') = hu 5(x,x') + I dydy' P Q (x, y) - 6 _ - x (B.2) 

{(EW)-«.)ft(y.^- 

qet/ 

(E^(yV q (y))Po(y',x')}. 

Here P Q is the projection operator in the first Landau level and the direct and exchange 
contributions are determined by the set U formed by the momenta associated to the selected 
filled states. At this point, let us assume that U = F. From the properties of the reflections 
it follows that 

R J H F HF \x, x') f(x') dx! = J H F HF) (Rx, Rx')R /(x') dx!. (B.3) 
Noting the invariance property that follows from Eq. A. 8, 

p (Rx,Rx') = J2MRy) vU R y') 
= E </Wy) v^ q (y') 

qGB 
= P (X,X') 

the invariance of the Hamiltonian follows, 

2 

H F HF) (Rx, Rx>) = hcu 5(x, x') + J dydy' P (Rx, Ry) — e _ — x (B.4) 
{( E ^(Ry'^Ry') - n ) P (Ry, Rx') - 

q£F 

(EvW^DW.*')} 

qeF 

= < F) (x,xO, 
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where in the sums again use of A. 8 has been made, as well as the fact that if q e F then its 
reflection about D is also in F. Using the above property in Eq. B.3 one gets 



R J #^ F) (x,x') /(x') dx! = J H { F HF \x,x')R /(x')dx' 

/From the arbitrariness of / then follows the commutativity of the reflection operator with 
H( HF \ 



[H ( F HF) ,R]=0. 



and consequently 



RHf^ | <p p > = e(p) R\<fip) 
= H F HF) R\<p p ) 
= e(-Rp) I Vrp ), 

leading to the symmetry of the spectrum with respect to reflection about the line D defined 
byd = (l,l), 

e(p) = e(Pp). (B.5) 

Clearly, the same conclusion arises with respect to reflections on the line D defined by 
d = (-l,l). 

The next task is to consider symmetry properties under a translation in the special vector 
a defined in Appendix A . In this sense, it is worth noting that if p e F then p + 5 E B — F 
and moreover, R p exactly gives p + 5 under a reflection about the line defined by the 
boundary of F and B — F. As we shall see, the relation between values of the self energies 
under such operations will allow us to show the property we are after. 

We start out by considering the states in the sector B — F and exploit particle hole 
symmetry. Adding the hamiltonian obtained through such states to that given in B.4 one 
obtains the Hartree-Fock form appropriate for a filled Brillouin zone, 



if <T(x, x') = H F HF \*, x') + #£?(x, x ) (B.6) 

/e 2 
dydy' P Q (x,y) - _ — 

£^(y'K(y) p (y',x') 



qS-B 

e 2 



hu 5(x, x') - J dydy' P (x, y)- 

P (y,y')Po(y',x'). 



y-y' 



It can be readily verified that the kernel H B commutes with all the magnetic translations, 
that is, 



[T a ,H ( » F) } = 0. (B.7) 



12 



Henceforth, since in particular it commutes with all the Tg leaving invariant the lattice L, it 

( IT JP\ 

follows that the ip p are eigenfunctions of H B . Additionally, they all have exactly the same 
eigenvalue. This can be verified by considering a translation in any vector a' = r 2 n x 8' 
commuting with it, which changes the momenta in an also arbitrary quantity 5' (Appendix 
A). One has H B HF ^ \ <p p ) = e B (p) \ ip p ),that may also be written as H B HF ^T r 2 nxS , \ 
VpS' ) — e B(p)T r % nX 5> | v?p-<5' )• Multiplying this latter expression from the left by the 
inverse translation one gets H B HF ^ \ ip p s' ) = €b(p) I Pp-S' ), which also equals e^(p — 5') \ 
ippS' )• Therefore, e#(p) = e_e(p — 5') = 6b, a constant. This property enforces the 
proportionality of H B with the projection operator P . At this point it is worthwhile 
to recall that since we are considering functions on the first Landau level, the Dirac delta 
function appearing in the definition of the HF kernels is equivalent to the projection operator. 
Then, 

#f F) (x,x') = hu 6(x,x!) + e L P (x,x') (B.8) 

Thus, all the basis functions are eigenfunctions of the Hartree-Fock hamiltonian associated 
with the filled Landau level H B HF ^ \ ip p ) = e B \ ip p ), with the same eigenvalue eg = Uu} +6l. 

The next point to consider is the connection between the kernels Hp 1 ^ and H B H ^p\ It 
follows from expressing the basis functions with momenta in B — F by means of translations 
of the basis states with momenta in F, and considering the transformation properties under 
the special translations T r 2 nxS of the projection operator and the Coulomb potential, as 
discussed in Appendix A. It follows 

#£?(x, x') = hw 8(x, x') + J dydy' P G (x, y) |y f y/| (B.9) 

{ E T ri nx^q(y') %l „ x5 ^(y')]* ~ U ) P a (y, x') - 

E T rl nx5 V?q(y) [^2 nx 5^q(y')]* Poty 
q€F 



hu>o5M + [ dydy' P fay)J%) nxS -J— 



y-y 

{( E ^(y')My') - n a ) nX(5 p ( y ,x') - 

qg-F 

.(y') 



E^q(y)^(yO ^ nx5 P (y',x')}, 

q€F 

which may be expressed in the more compact Dirac notation as, 

H B H _ F J = T r 2 nX(5 H F HF ^T_ r 2 nx5 . (B.10) 

Physically, this relation means that the HF state associated to F is simply the space trans- 
lation of the state associated to B — F. Therefore, within the first Landau level the following 
relation holds, 

£b = H F HF ^ + T r 2 nx 5 Hj^^T-j.* nx<5 , (B-H) 
which after acting over a basis function produces 
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e B = e(p) + e(p-5), 



(B.12) 



The sum of the energies at p and the shifted value is strictly constant. Using the symmetry 
of the spectrum under reflection it also follows that 



Thus, the sum of the energies associated to the states which have momenta related by a 
reflection in the boundary of the regions F and B — F is also exactly constant. Finally by 
selecting p = Pf , with p^ in the boundary between both regions, and noticing that in that 
case Pf = R(pf — 8), it follows that, due to the continuity of the spectrum (absence of gaps) 
the energy at any point of the boundary of the filled states strictly equals e(pp) = ep/2. 
Therefore, the criterium for a minimum given in reference" is satisfied, and the Hartree-Fock 
solution with all states in F has a local minimum of the energy. 



e B = e(p) + e(R(p-5)). 



(B.13) 
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